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ABSTRACT 

We present an analysis of the polarization of compact radio sources from six pointings of the West- 
erbork Synthesis Radio Telescope (WSRT) at 350 MHz with 35% coverage in A 2 . After correcting 
for the off-axis instrumental polarization with a simple analytical model, only a small number of 585 
strong sources have significant polarizations at these wavelengths. The median depolarization ratio 
from 1.4 GHz for the strongest sources is <0.2, reinforcing the likelihood that radio galaxies are found 
in magnetized environments, even outside of rich clusters. Seven sources with significant 350 MHz 
polarization were selected for a more in-depth Faraday structure analysis. We fit the observed values 
Q/I and U/I as a function of A 2 using both a depolarizing screen and two component models. We 
also performed RM Synthesis/Clean and standard fitting of polarization angle vs. A 2 . We find that 
a single rotation measure (RM), as found using polarization angle fitting or simple screen models, 
commonly provides a poor fit when the solutions are translated back into Q, U space. Thus, although 
a single "characteristic" rotation measure may be found using these techniques, the Faraday structure 
of the source may not be adequately represented. We also demonstrate that RM Synthesis may yield 
an erroneous Faraday structure in the presence of multiple, interfering RM components, even when 
cleaning of the Faraday spectrum is performed. We briefly explore the conditions under which rota- 
tion measures and Faraday structure results can be reliable. Many measurements in the literature do 
not meet these criteria; we discuss how these influence the resulting scientific conclusions and offer a 
prescription for obtaining reliable RMs. 

Subject headings: Galaxies: intergalactic medium - galaxies: clusters: intracluster medium - polar- 
ization - radio continuum: galaxies - techniques: polarimetric 



1. INTRODUCTION 

By characterizing the Faraday structure in radio 
synchrotron sources, properties of the magneto-ionic 
medium can be probed, such as magnetic field strength 
and orientation, as well as distribution of the relativistic 
and thermal electron populations. Radio arrays such as 
the Westerbork Synthesis Radio Telescope (WSRT), the 
Expanded Very Large Array (EVL A) , the Low Frequency 
Array (LOFAR), the Australia Telescope Compact Ar- 
ray (ATCA), the Allen Telescope Array (ATA), and the 
planned Australian Square Kilometer Array Pathfinder 
(ASKAP) are well suited for Faraday structure studies 
due to their enhanced A 2 sampling capabilities, e.g., wide 
relative bandwidth (AA 2 /A 2 rai „), and high spectral reso- 
lution (c5A 2 ). 

We represent the complex linear polarization by 



ties. We use q, u to represent the fractional values Q/I, 

Traditionally, most polarization studies have deter- 
mined rotation measures by fitting 



X(A 2 ) = xo + A 2 RM, 



(4) 



P = Ipe 2lx = Q + iU = I(q + iu) 



(1) 



where p and x are the degree and angle of polarization, 
given by 
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(3) 



and I, Q, U are the Stokes parameters for the total and 
orthogonal components of the linearly polarized intensi- 



where RM is the Faraday rotation measure, with little 
or no attention paid to the behavior of the fractional po- 
larization. A common practice has been to restrict RM 
fitting to regions of A 2 space where p(X 2 ) is constant or 
decre ases monotonically (e.g., iSimard-Normandin et al.l 
1981), which would occur for a foreground rotating or 
depolarizing screen. This is sometimes done even when 
data showing a rise in p(X 2 ) at shorter wavelengths exists, 
ignoring evidence that multiple RM components may be 
present. Others restrict their fitting to A < Ai / 2 (de fined 
by p(Ai/2)/p(0) = 0.5), beyond which iBurnl (|1966[ ) sug- 
gests that si gnificant no n-linear behavior in x(A 2 ) i s cx ~ 
pected (e.g., [Haves 1975?). In other cases, significant non- 
linea r behavior in y(A 2 ) is observed (e.g.. lMorris fc Bergd 
11964 IRqv et all 12005). but no modeling of this anoma- 
lous behavior is made and the RM from the poor linear 
fit to x(A 2 ) is reported. Others require that p is above 
some threshold and/or that a minimum signal-to- noise 
value is present in the observations but do not report the 
behavior of p(A 2 ), which may hold inform ation regard- 
ing the underlying Fara day structure (e.g.. [Clarke et aLl 
2001, Br own et al.ll2007l ). 
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The only situation where dx/dA 2 and p(X 2 ) are con- 
stant, and therefore unimportant in determining the 
Faraday structure, is when there is a single uniform Fara- 
day screen completely in the foreground. In all other 
cases, including all or most physically realistic ones, 
more sophisticated m odeling is required. For example, 
iFletcher et al.l ([2004D consider both polarization degree 
and angle i n their study of the mag netic field o f M31 . 
In addition. iRossetti et al.l (|2008t ) and lFanti etall (|2004f ) 
employed simple models of depolarization and x rota- 
tion to examine Compact Steep-Spectrum (CSS) sources 
at A < 21 cm. We will briefly summarize some of the 
classic models where p{\ 2 ) 7^ constant. For a detailed 
di scussion on depolariz ation effects, we refer the reader 
to ISokoloff eFalT (fl998h . 

If the thermal electrons are spatially coincident with 
the relativistic, synchrotron emitting electrons, e.g., then 

, «. sin(A 2 F c ) , , 

P(A 2 ) oc -^f>. (5) 

as in the uniform slab model of IGardner fe Whiteoakl 
(1966), where F c is the Faraday depth th rough the slab; 
F c can be thought of as an "internal" RM. ICioffi fc Jonesl 
( 1980) showed that the observed depolarization and \ ro- 
tation can have considerable differences depending on the 
geometry assumed, even for simple cases such as cylin- 
ders and spheres. 

For a foreground screen consisting of many unre solved 
compo nents with a random distribution of RMs, Burn 
(1966) modeled the observed fractional polarization as 

p(X 2 ) ex exp(-2a 2 , M A 4 ) (6) 

where o\ M is the variance of a Gaussian dispersion in 
RM across this so-called "mottled" screen. Modifications 
to this model have been proposed, e.g., by IRossetti et al.l 
(2008) who include the effect of filling factors. 

Two interfering foreground RM components will also 
produce non-A 2 behavior in angle, and changes in frac- 
tional polari zation that can ris e or fa ll with increasing 
wavelength. [Goldstein fc Reedl (|1984[ ) describe the ob- 
served polarized flux from two such components as: 

Pobs = Pi[l + k 2 + 2fccos( X i - X2)] 1/2 (7) 

where k = P2/P1 < 1 is the ratio of the the polarized 
fluxes, xi an d X2 are the polarization angles at the ob- 
servation frequency. 

With adequately sampled data in A 2 space, all of the 
above cases can, in principle, be distinguished. In prac- 
tice, however, A 2 sampling is inadequate to map out the 
Faraday structure, and even the large fractional band- 
widths of the WSRT or the EVLA can be insufficient. 
As we will illustrate below, determination of the Faraday 
structure requires observations which detect the varia- 
tions in both p{\ 2 ) and dx/dA 2 . In particular, the result 
dx/dA 2 « constant can occur over a substantial range in 
A 2 even with underlying Faraday structure. Whether or 
not failure to diagnose the presence of underlying Fara- 
day structure is acceptable depends on the particular sci- 
entific goals, as we discuss further below. 

In Section [5] we present our WSRT observations at 350 
MHz and the determination and removal of the off-axis 
instrumental polarization. We present the results of our 



polarization and Faraday structure analyses in Section 
[3l There we characterize the depolarization of our sam- 
ple of 585 compact sources and give a brief overview of 
the polarization diagnostics and Faraday structure mod- 
eling employed. We then detail the modeling results on 
seven sources with significant 350 MHz polarization and 
the discrepancies between fitting q{\ 2 ) and m(A 2 ) and 
other other techniques. In Section|4]we use the results of 
some simple experiments to demonstrate some of the in- 
adequacies of common RM determination methods such 
as x(^ 2 ) fitting and RM Synthesis. We also offer some 
recommendations for reliable RM determinations. A dis- 
cussion of our findings, including the science implications 
of RM ambiguities, is presented in Section O 

2. WESTERBORK 350 MHZ OBSERVATIONS AND 
INSTRUMENTAL POLARIZATION 

2.1. Observations and Data Reduction 

We observed six fields with the WSRT in 2008 and 
2009, originally selected for possible large-scale diffuse 
polarization fou nd in the NRAO Ver y Large Array Sky 
Survey fN VSS. iCondon et alJ IT9981 through a repro- 
cessing by iRudnick fc Brownl ([2009) . To minimize the 
contribution of polarized Galactic foreground emission, 
we have selected fields with |6|>42°. Observations were 
made in spectral line mode with a central frequency of 
345 MHz, 70 MHz bandwidth, and 1024 channels over 
8 intermediate frequency (IF) sub-bands, yielding full 
Stokes parameters. Even though the central frequency 
is 345 MHz, we will continue to refer to this band as 
the 350 MHz band to comply with the established con- 
vention. Various array configurations were used and are 
shown in Table [TJ The nominal synthesized beam size 
varies with array configuration and observing frequency, 
but is approximately 70" in RA for our observations. 
Due to the East- West array configuration, the beam be- 
comes elongated in the North-South direction by a factor 
of csc(<5). 

We will now summarize the key elements of the data 
reduction and calibration process; for a complete des crip- 
tion we refer the reader to lBrown fc Rudn ick (2009J). All 
reduction was done using standard techniques in AIPS, 
correctly accounting for the WSRT linearly polarized 
feeds, and including several iterations of amplitude and 
phase self calibration for total intensity. Flux calibra- 
tors were observed, and are listed in Table [TJ The AIPS 
procedure LPCAL was used to correct for polarization 
leakage between the X and Y orthogonal linear polar- 
ization receivers. Additionally, calculation of a time- 
independent phase correction between linear polariza- 
tions X and Y was attempted for each channel using a 
polarized calibrator observed during the run. The polar- 
ized calibrator 3C345 was used for all fields except Field 
B, for which DA240 was used instead. Unfortunately, 
a solution was not found for every channel, rendering 
those channels without a solution useless for Stokes U 
measurement. 

Cleaning and imaging were also done in AIPS, where 
4°x4° images in Stokes /, Q, and U were created for each 
channel. The community is just beginning to experi- 
ment with the much simpler problem of multifrequency 
synthesis/cleaning in total intensity, where one or two 
spectral parameters can be used to characterize the fre- 
quency dependence, and the biases there have not yet 
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TABLE 1 

Summary Of WSRT 350 MHz Observations 



Field 


RA 


DEC 


6 1 


Array 


Exposure 


Common 


Calibrators 




(J2000) 


(J2000) 


C) 


Configuration 


(Hours) 


Beam (") 




Aries-Pisces 


01:09:14.30 


+13:09:58.0 


-49 


Mini-short 


12 


325 x 70 


3C147, 3C295 


Coma SW 


12:54:08.00 


+26:42:00.0 


+89 


Special 2 


12 


70 x 70 


3C147, 3C295 


Coma NW 


12:54:08.00 


+27:58:00.0 


+88 


Special 2 


12 


70 x 70 


3C147, 3C295 


Coma NE 


12:59:52.00 


+27:58:00.0 


+88 


Special 2 


12 


70 x 70 


3C147, 3C295 


A 


14:53:00.00 


+40:25:00.0 


+62 


Maxi-short 


1.2 


125 x 70 


3C147, 3C48 


B 


16:20:00.00 


+60:12:00.0 


+42 


Maxi-short 


12 


105 x 70 


3C48, 3C295 



Approximate Galactic latitude at field center. 
2 Special array configuration is 36m+54m+72m+90m. 



TABLE 2 

Number of 78 kHz Channels Used For IF Band Averaging 



Field 


Stokes 


IF1 


IF2 


IF3 


IF4 


IF5 


IF6 


IF7 


IF8 






(376.4) 1 


(367.7) 


(358.9) 


(350.2) 


(341.4) 


(332.7) 


(323.9) 


(315.2) 




I 


100 


97 


68 


79 


100 


99 


64 


32 


Aries-Pisces 


Q 


100 


99 


76 


79 


100 


99 


64 


32 




u 


69 


96 


58 


24 


85 


28 


64 


32 




i 


101 


101 


59 


62 


101 


98 


96 


92 


Coma SW 


Q 


101 


100 


58 


69 


100 


97 


97 


97 




u 


86 


52 


[7] 2 


40 


96 


38 


39 


59 




I 


101 


101 


70 


96 


101 


99 


99 


93 


Coma NW 


Q 


101 


100 


73 


79 


99 


97 


98 


96 




u 


91 


84 


21 


75 


100 


27 


78 


32 




I 


101 


101 


56 


29 


101 


99 


91 


87 


Coma NE 


Q 


101 


100 


64 


62 


100 


98 


96 


95 




u 


98 


85 


[5] 2 


16 


99 


54 


68 


57 




i 


92 


101 


71 


96 


96 


99 


101 


88 


A 


Q 


95 


101 


76 


96 


96 


99 


101 


88 




u 


30 


[I] 2 


26 


78 


65 


[2] 2 


17 


88 




i 


86 


94 


42 


41 


100 


92 


90 


70 


B 


Q 


86 


97 


64 


76 


101 


97 


101 


93 




u 


52 


98 


76 


88 


101 


97 


101 


69 



1 IF band central frequencies are g 
2 Bracketed values identify IFs where U averagin; 

been characterized. Q and U have much more complex 
behavior as a function of frequency and will require ex- 
tensive experimentation in the future. Therefore, each 
channel and Stokes quantity was cleaned separately in 
AIPS with IMAGR using a loop gain of 0.1 and 15,000 
clean components per field. Images of Stokes V (circular 
polarization) were made to verify that no leakage into V 
was present, under the assumption that it is negligible for 
typical astrophysical sources. Typical channel sky RMS 
values of ^3-5 mJy/beam were obtained in the cleaned 
/ images, and ~l-3 mJy/beam for the Q, U images. All 
images for a field with sky RMS < 5 mJy/beam (uncor- 
rected for primary beam attenuation) were convolved to 
a common beamsize, allowing channel averaging to be 
performed as described below. 

Average images of Stokes I, Q, and U for each of the 
eight IFs were constructed from the individual channel 
images, along with a total intensity map averaged over 
all eight IFs. The number of channels used for each band 
average image is listed in Table [2J Channels with imag- 
ing problems, such as strong artifacts due to radio fre- 
quency interference (RFI), were excluded. In addition, 
U imaging was not performed on channels where no X- 
Y phase correction was found. We supplemented these 
data using images from the NVSS to provide measure- 
ments of Stokes /, 0, U at 1.4 GHz. The NVSS images 



iven below IF number in MHz. 

; was not performed due to too few channels. 

were convolved to the corresponding WSRT field's beam- 
size. See Table [T] for an overview of the field properties, 
including common beam convolution sizes. The Coma 
fields were imaged using a (somewhat smaller than nom- 
inal) restoring beam of 7 0"x70" as part of another study 
(jBrown fc Rudnicldl20Tl . 

Total intensity images of the six fields are shown in 
Figure [TJ Images of the linear polarization at RM=0, 
taken from the results of RM Synthesis (see Section 

I3.2.1j) . are shown in Figure [2] Note that the polarization 
maps at RM=0 are pervaded by diffuse Galactic emis- 
sion ( e.g.. iBrentiens fc de Bruvn 2005. Schnitzcle r et al.l 

[200l IWolleben et al.ll2010l iBernardi et al.ll2010D ~ 

2.2. Instrumental Polarization of WSRT at 350 MHz 

To identify sources with either real or instrumental po- 
larization, we first selected sources in each field with I/cri 
> 30 in the all-IF Stokes / image, yielding 585 total 
sources for the six fields. We then extracted Stokes /, Q, 
U from each of the WSRT individual IF and NVSS im- 
ages, at the peak location in the all-IF Stokes I image. A 
background subtraction was performed for each measure- 
ment using a rectangular region about the source, of in- 
ner dimension lx the synthesized beam dimensions and 
outer dimension 2x the synthesized beam dimensions. 
The RMS deviation within each annulus was adopted as 
the statistical error in each measurement. 
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Fig. 1. — Total intensity images of the six fields observed at 350 MHz with the WSRT. Top row, from left: Aries-Pisces, Field A, Field 
B. Bottom row, from left: Coma SW, Coma NW, Coma NE (See Table[T}. Images are 4x4 degrees. Diffuse emission from the Coma halo 
and relic are visible in the Coma images. Also visible are residual imaging artifacts near the strongest sources, common for the WSRT. 
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Fig. 2. — Linear polarization images at RM=0 of the six fields, taken from the results of RM Synthesis, observed at 350 MHz with the 
WSRT. Top row, from left: Aries-Pisces, Field A, Field B. Bottom row, from left: Coma SW, Coma NW, Coma NE (See TableQ}. Images 
are 4x4 degrees. Note the diffuse Galactic emission which pervades each field. Also visible are residual imaging artifacts near the strongest 
sources, common for the WSRT. 
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TABLE 3 

Model Fit Parameters for Instrumental Polarization 



1000 2000 3000 4000 
R (arcsec) 



5000 



6000 



Fig. 3. — Plot of average p cor r at 350 MHz for sources with at 
least six of the eight IF measurements satisfying Pmeas/o~p > 1.2 
before instrumental correction has been applied. The instrumen- 
tal polarization increases with off-axis radius. The vertical line at 
4500" corresponds to the radial limit of our instrumental polariza- 
tion model fitting. 

For the purpose of illustrating the instrumental po- 
larization, we apply the simplest bias correction to the 
polarization amplitude: 



P2 _ a 2 

1 meas " Pi 



(8) 



which i s an approximation to the " most probable estima- 
tor" of iWardle fc Kronbergl (|1974T >. good for P corr /r7p > 
0.5. This most prob able estimator is the best available 
for Pcorr/vp > 0.7 (|Simmons fc Stewartlll985|) . and we 
only report results well above this limit. We use a prop- 
agated error calculation for up based on the observed 
errors in Q and U . 

Figure [3] shows p CO rr = Pcorr/ I, averaged over mul- 
tiple IFs, vs. off-axis radius for the WSRT data set, 
illustrating the instrumental enhancement of fractional 
polarization with rad i us as mentioned previously by 
Ide Bruvn fc Brentiensl (120051) an d investigated at 1.4 
GHz bv IPopping fc Bra un (2008). In this work, we de- 
termined the instrumental polarization behavior in both 
Q and U for each IF in order to perform a first order 
correction. For each IF, we first selected from the 585 
initial sources those satisfying P meas /crp > 2 and plot 
q and u as a function of their locations relative to the 
pointing center (Figure |4]). In several IFs there is a clear 
quadrupole pattern, in general possessing a greater mag- 
nitude in q than u. The observed q quadrupole pattern 
is oriented coincident with the orientation of the X and 
Y linear dipole feeds on the WSRT, which face the sky 
perpendicular to each other and form Stokes Q by the 
linear combination of XX* and YY* . The observed u 
quadrupole pattern, which is formed from a linear com- 
bination of XY* and YX*, is offset 45° on the sky with 
respect to the q pattern, as one would expect. 

To quantitatively model the instrumental polarization 
for each IF, we made the following cuts to the data. 
Outside 4500" the instrumental polarization rises sharply 
and we do not attempt any correction beyond that limit, 
cutting the total number of sources from 585 to 3353. 
We then required P m eas/vp > 2, yielding roughly 100 
sources per IF. For each IF we fit a double cosine func- 

1 This includes 36 duplicated sources observed at different off- 
axis positions due to the multiple pointings for the Coma field. 
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fxlO -3 /") 


(°) 


1 


O 
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47 


5 


Q 
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90 




u 
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33 


6 


Q 
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91 




u 
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0.7 


4 


7 


Q 


0.15 


0.78 


85 




u 


0.13 


0.13 


16 


8 


Q 


0.0005 


1.9 


108 




u 


0.06 


0.65 


62 



tion to each set of q and u of the form 
f(r, PA) = Ae Br cos(2PA H 



C), 



(9) 



which includes the distance from the pointing center, 
r, and position angle, PA, of the source. This yielded 
16 total sets of parameters, which are given in Table |3l 
We then produced corrected Q, U observations for each 
source by subtracting the modeled instrumental contri- 
bution. 

The instrumental polarization is weak near the point- 
ing axis, generally much less than 1%, but grows to as 
much as 6% in q near R p b (half power radius of the pri- 
mary beam) for the odd numbered IFs. In u the instru- 
mental contribution is <3% at R p b for all IFs. 

By examining Figure [5j one can see evidence o f the 1 7 
MHz modulation, as found by IPopping fc Braun] ((2008) , 
in the q models for IFs 1, 3, 5, and 7, which are separated 
by ssl7 MHz. In these IFs, the instrumental polariza- 
tion is stronger by factor of roughly 2-3 at R p b over the 
neighboring even numbered IFs. This effect is much less 
pronounced in u, as seen in Figure [5j 

After correction for instrumental polarization, there is 
still a significant polarized flux bias from a variety of 
factors which differ from one IF to another, including 
the noise bias (including random noise and residual side- 
lobe structures) and non-quadrupole components to the 
instrumental polarization as a function of IF and two- 
dimensional location within the primary beam. These 
are not well modeled by Equation [SJ so in order to make 
a practical model for the polarized flux bias we took an 
empirical approach and measured the median (p me d35o) 
and RMS scatter (p SC atter35o) among the 350 MHz IFs 
of the polarized fraction for each of 335 sources with 
r < 4500". We expect that residual instrumental po- 
larizations, sidelobe structures, and noise will all vary 
from IF to IF, and that p sca tter350 will therefore provide 
an estimate of all of these contributions. On the other 
hand, p m ed350 provides an estimate of the true polarized 
flux, along with a bias related to p SC atter350- These are 
plotted vs. each other in Figure [SJ Different symbols 
represent different levels of NVSS polarized flux for the 
same sources. 

There was no significant difference in the distribution 
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Fig. 4. — Plot of q and u measurements for sources with P bs/<^P > 1.2 from all six fields as projected on the sky, illustrating the 
radial and azimuthal behavior of off-axis instrumental polarization. Blue points are negative, red points are non-negative; the point size is 
related to the magnitude of the measurement. The quadrupole pattern described in the text is apparent. The black circle, of radius 4500", 
surrounds the region used for model fitting of instrumental polarization. 



as a whole between sources with no NVSS polarization 
and sources with moderate or strong NVSS polarization. 
Therefore, the bulk of p m ed350 values are likely due to the 
instrumental contributions described above, as opposed 
to intrinsic polarizations. We fit the distribution and 
found 

Pmed350 = Pscatter350 

) 2 + (0.006) 2 (10) 

We then adopted this calculated value as the effective 
polarization bias to be subtracted in quadrature from 
each of the measurements when doing statistical anal- 
yses. If a source had an intrinsic polarization equal to 
1.5xp scat ter350 which would add in quadrature to the cal- 



culated value of Pmed350, then the source would be found 
on average at the dotted line in Figure [5] Only three 
sources out of 335 exceed this value (and at least two do 
have well-behaved polarization behaviors), so we adopt 
this as our upper limit for the purposes of calculating 
depolarization ratios. 

We note that changes to the empirical fit in Equation 
[TU]will have a small effect on the statistical analyses in 
which it is employed. For example, if the fit value of 
Pmed350 is overestimated the above procedure may elim- 
inate some sources that have significant real structure in 
p(X ). However, the number of such sources is small, as 
discussed further below, so we ignore that issue in order 
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Fig. 5. — Model fits to the instrumental polarization of WSRT at 350 MHz as projected onto the sky (North is up, East is left). Top 
row: q for IF1 - IF4. Second row: q for IF5 - IF8. Third row: u for IF1 - IF4. Bottom Row: u for IF5 - IF8. The frequency dependence 
can be seen in q by noting the increased amplitude in the odd numbered IFs. The black circle is of radius 4500". 



to examine the depolarization properties of the sample 
as a whole in the following section. Since the residual 
bias correction from Equation [10] is not applied in the 
individual source modeling described in Section [3.2.31 it 
has no effect on the outcome of those analyses. 

3. RESULTS 

3.1. Polarization Properties of the Overall Sample 

Starting with the sample of 335 sources discussed 
above, we determined their polarized fluxes in the 1.4 
GHz NVSS survey. We first convolved the NVSS /, Q 



and U images to the same beamsize as used in each cor- 
responding WSRT field, then measured the I, Q and U 
fluxes at the locations of each total intensity peak in / at 
350 MHz. Background subtraction and error estimation 
were performed using the same rectangular region about 
the source as described previously for our WSRT mea- 
surements. We then calculated the polarized flux (and 
fractional polarization) after correcting for the noise bias, 
according to Equation [H After correction, we found that 
102 of the 335 sources had significant polarizations at 
1.4 GHz {PcorriA/o-piA > 2), and for each we calculated 
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Fig. 6. — Plot of the median of p in the eight WSRT IFs vs. 
the RMS scatter of p among the IFs, to determine the polariza- 
tion bias remaining after correction for instrumental polarization. 
Circles represent sources exhibiting no polarization in the NVSS. 
Xs represent sources exhibiting moderate or strong polarization in 
the NVSS. Solid triangles show sources chosen for modeling as de- 
scribed in the text. The solid line is the best fit for all bands for 
each of the 335 sources used to model the instrumental polariza- 
tion, and the dashed line shows the defined upper limit discussed 
in the text. 
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Fig. 7. — Plot of bias corrected median P350 vs. bias corrected 
pi.4 for 102 sources with significant polarization at 1.4 GHz, de- 
fined as P/crp > 2. Xs represent upper limits at 350 MHz for the 
median fractional polarization. Solid symbols represent the sources 
we modeled, except for NVSS J162740+603900 which did not have 
a significant detection in polarization at 1.4 GHz. Circles represent 
sources whose median 350 MHz values are formally upper limits, 
although they were clearly detected in some IF bands. Diamonds 
are significant detections at both bands, shown with their errors. 

the upper limit to their polarized fractions at 350 MHz. 
These are plotted in Figure [3 

The upper limits on the 350 MHz polarized fractions 
are largely independent of the fractional polarizations at 
1.4 GHz. There is a rough upper limit to the distribution 
visible in Figure [7] likely due to the fact that at low fluxes, 
only large values of pi.4 can be detected, and the upper 
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Fig. 8. — Plot of the upper limits to depolarization ratio from 
1.4 GHz t o 35 MHz vs. pi.4 for the set of 335 sources described 
in Section l3.ll Depolarization ratio is defined as P35o/pi.4- 

limits on P350 will therefore also be high. Lines of unity 
slope on this diagram indicate specific depolarization ra- 
tios (j>35o/pia)- Upper limits to the depolarization ratios 
vary from <0.03 to <2, with a median of <0.3 . 

In Figure [8] we plot the median upper limit to the de- 
polarization ratio as a function of P1.4. The decreasing 
values indicate the observational bias that we can only 
measure low depolarization upper limits for the highest 
values of pi. 4. The median upper limit for the 20 highest 
P1.4 sources is ~0.2. A more conservative requirement of 
PcorriA/ '&pia > 3 did not change the overall distribution 
of the depolarization ratios. 

3.2. Model Fitting of Individual Sources 
3.2.1. Source Selection 

From the set of 335 sources (with r < 4500" and 
I/crj > 30 at 350 MHz, and disregarding pia), we 
identified a subset based on their Faraday Dispersion 
Function (FDF) using the RM Synthesis technique 
(jBrentiens fc de Bruynll2005h . This allows for the best 
signal to noise averaging of all the data, since Q{\ 2 ) and 
U (A 2 ) can be summed as vectors after correcting for each 
assumed RM. 

The obser ved FDF, F ((/)), is construc ted (using the 
formalism of iBrentiens fc de Bruvnl 12005) thusly: 

N 

Fit) = F{4>) * R{<t>) =KY, wtPie-W*-*) (11) 

i 

N 

R(cj>)=K y £ / w i e- 2i ^-^ (12) 

i 

K =\H w A ( 13 ) 



at an arbitrary Faraday depth, <f>, which replaces the 
usual rotation measure; in practice one chooses a range of 
Faraday depths to reconstruct a Faraday spectrum. The 
quantities Pi and Wi are the observed vector polarization 
and applied weight, respectively, at locations of sampled 
A 2 . The quantity A 2 , is the mean A 2 of the set of ob- 
servations, and the reconstructed FDF is represented at 
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TABLE 4 

RM Synthesis Capabilities of WSRT 



Frequency 


<A 2 > 


AA 2 


A 2 

min 


«5A 2 


S<j> 


max-scalc \4>max\ 


(MHz) 


(m 2 ) 


(m 2 ) 


(m 2 ) 


(m 2 ) 


(rad/m 2 ) 


(rad/m 2 ) (rad/m 2 ) 


310-380 


0.76 


0.31 


0.62 


3.1 x 10~ 4 


12 


5.0 5700 



TABLE 5 
Table of Modeled Sources 



Source 


RA 


DEC 


Off-axis 


Pos. Angle 




(J2000) 


(J2000) 


Radius (") 


(°) 


NVSS J010616+125116 1 


01:06:16.8 


+ 12:53:22 


2786 


250 


3C33S 2 


01:08:50.7 


+ 13:18:43 


649 


326 


NVSS J011136+132437 


01:11:36.2 


+13:25:41 


2268 


65 


NVSS J011204+124118 


01:12:04.5 


+12:42:39 


2962 


123 


NVSS J125630+270108 


12:56:30.5 


+27:01:10 


3816 


150 


NVSS J162408+605400 


16:24:08.8 


+60:54:04 


3134 


35 


NVSS J162740+603900 


16:27:41.0 


+60:39:05 


3783 


63 



Resolved as double source in unconvolved NVSS image (with NVSS J010615+125210) 

2 NVSS J010850+131831 



A 2 = Aq. Note that the actual F(<f) is obtained by decon- 
volving the Rotation Measure Spread Function (RMSF, 
R((f>)), which is the normalized response function in Fara- 
day space, from the observed F(cj>). We briefly discuss 
the deconvolution procedure, RM Clean, in Section l3.2.3l 
Unless otherwise noted, all FDFs and RMSFs in this 
study were constructed using uniform weighting. For this 
paper we use an over-tilde to represent transformed po- 
larization quantities unless otherwise noted, e.g., P rep- 
resents the magnitude of the FDF, Q represents the real 
part of the FDF, and so on. 

We used all channels where sky noise in Stokes Q and U 
were < 5 mJy/beam (uncorrected for primary beam at- 
tenuation), with the number of channels listed in TableO 
A typical RMSF for the WSRT 350 MHz band is shown 
in Figure [101 this RMSF was constructed using roughly 
400 channels across the full band. The main lobe of each 
RMSF had a characteristic FWHM -12 rad/m 2 . Nom- 
inal RM Synthesis capabilities of the WSRT 350 MHz 
band are given in Table |4] We note that no instrumental 
polarization correction has been applied to the data used 
to construct these FDFs, since the corrections were de- 
termined only when the channels were averaged within 
each IF band. 

A coarse initial search over Faraday depths between 
±1000 rad/m 2 was performed at a resolution of 10 
rad/m 2 . Once we had determined that no significant 
power existed outside a Faraday depth of ±200 rad/m 2 , 
we performed a finer search between ±200 rad/m 2 with 
a resolution of 1 rad/m 2 . 

To make an initial cut to the set of sources, the loca- 
tion {4>max) and amplitude (P ma x) of the peak in P(<ft) 
were determined for each source, and those with a peak 
amplitude of P ma x > 3 mJy/beam/RMSF (uncorrected 
for primary beam attenuation) were selected. In all, 116 
of the original 335 sources passed this criterion, with a 
minimum signal to noise in P max of 3.4. This removed 
many sources from the sample whose observed polariza- 
tion may be enhanced artificially, e.g., by noise, which 
places power at all Faraday depths in the FDF. We note 
that many of these remaining detections are due to in- 
strumental polarization which is not corrected in the all- 



channel FDF. 

For each of these 116 sources, we then examined the 
IF averaged Q and U measurements, corrected for instru- 
mental polarization. Because sources could have different 
fractional polarizations for different IFs, we did not de- 
mand that they have strong signals in all IFs. Sources 
with at least four of the eight IFs satisfying Pjap > 2, 
and U/ajj > 4 were then selected from the list of 116. 
We used only the U data for this cut because of the 
greater uncertainty in the instrumental correction for Q 
and the presence of occasional spuriously high Q values. 
Sources which exhibited a regular modulation in p(X 2 ) 
corresponding to to the —17 MHz modulation found by 
iPopping fc Braunl pOOl were excluded. All such sources 
were found beyond R ~ 4000", evidence of residual in- 
strumental polarization not fully accounted for by our 
model. Only three sources met all of these criteria. To 
those, we added four additional sources for modeling 
based on their high ratios of p me d35o/Pscatter350, putting 
them at or above the upper limit line shown in Figure 
[BJ These seven sources selected for modeling are listed 
in Table [5] along with selected properties from the lit- 
erature. Plots of q(X 2 ) and w(A 2 ) are shown in Figure 

m 

3.2.2. A Note on Bandwidth Depolarization 

The NVSS data at 1.4 GHz were constructed from 
two 42 MHz wide bands, centered at 1364.9 MHz and 
1435.1 MHz. Bandwidth depolarization for sources with 
|RM| < 50 rad/m 2 would yield (p -p)/p < 2% in the 
NVSS lCondon et al.lH998l) . Any source with |RMj high 
enough to suffer significant bandwidth depolarization in 
the NVSS would be severely depolarized in the IF aver- 
aged 350 MHz data, and would not have been selected 
for further investigation. 

3.2.3. Modeling Techniques 

As can be seen in Figure [9] the sources found to have 
sufficient polarization for modeling all showed structure 
in q(X 2 ) and w(A 2 ) inconsistent with a simple Faraday 
screen. A simple screen would result in sine and cosine 
waves in q and u with matched frequencies and ampli- 
tudes. In order to measure the Faraday structure of these 
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Fig. 9. — Observed NVSS + WSRT IF band averaged q(X 2 ) and w(A 2 ) for the seven sources modeled. Background subtraction and 
removal of WSRT instrumental polarization has been performed for each source, as described in the text. 
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sources, we therefore explored a variety of techniques. 
In particular, we used: A. Linear fit to x(A 2 ); B - R M 
Synthesis/Clean; C. Model fitting to q, u vs. A 2 using 
a two component foreground screen; and D. Model fit- 
ting using a single foreground screen with a mean RM 
and a separate depolarizing function. We have omitted 
an internal Faraday dispersion model for the following 
reason: internal depolarization in the Milky Way and 
nearby galaxies arises becau se the synchrotron a nd ther- 
mal plasmas are well mixed (Sokoloff et al. 1998). This is 
not true for extragalactic sources, where the depolariza- 
tion almost al ways arises wit h Faraday variations across 
the beam (e.g.. lTribblei ri991 and references therein). We 
now briefly discuss each of the models employed followed 
by the results. 

A. x(A 2 ). We determined the RM for each source using 
the most common method, minimizing the sum of the 
weighted residuals (i.e. chi-squared statistic, x 2 ) from 
fitting Equation 0] to the observed polarization angles 
x(A 2 ). There were often a number of different solutions 
with comparable values of Xmin based on our choices for 
the n7r ambiguities. We therefore made these choices to 
most closely match the results for the RM for the fore- 
ground depolarizing screen model described below. We 
calculated the errors in RM by the standard propagation 
of errors from the residuals to the fit, not from the errors 
in the original data points. Note that the reduced x 2 val- 
ues, xt = X 2 /dof (where dof — degrees of freedom), as 
listed in Tables I6l-fl2l are generally quite high, suggesting 
that these are not good fits, despite the apparently small 
derived errors in RM. 

B. RM Synthesis/Clean. For each of the seven sources 
the FDF was constructed using the q,u data, this time 
using the instrumental polarization corrected IF samples 
from the WSRT observations plus the NVSS data point. 
Uniform weighting for all A 2 samples was applied; we 
experimented with various weighting of the WSRT and 
NVSS samples used as input for RM Synthesis, but found 
negligible differences in the RM Clean solutions. A rep- 
resentative RMSF is shown in Figure [Til displaying the 
sidelobe structure due to the sparse A 2 sampling. The 
range of Faraday depth for the constructed FDF was this 
time limited to ±50 rad/m 2 , reflecting the maximum RM 
due to the A 2 separation of the IF averaged samples. 
The full-channel FDFs were first searched for significant 
power beyond ±50 rad/m 2 to ensure that this range of 
Faraday depth s was large enough. Our cu s tom version 
of RM Clean (|Brentiens fc de Bruvnl [20051 iHeald et al.l 
2009) was used to deconvolve the complex RMSF from 
the FDF, drastically reducing sidelobes and producing a 
more lucid representation of the Faraday structure. We 
used a gain factor of 0.1 and stopping criteria of either 
200 iterations or a peak to RMS ratio of 1.5 in the residu- 
als of p. These convergence criteria were found to strike 
the optimal balance between minimizing the residuals 
and limiting spurious clean components. 

The location (RM) and amplitude (po) of the two most 
dominant features in each cleaned FDF were extracted 
by Gaussian fitting to the cleaned Faraday spectrum. To 
reduce polarization bias, which would enhance the am- 
plitude of po solutions, we subtracted the mean of the 
residuals in p before performing the Gaussian fitting. As 
determined in some of our experiments, and also noted by 
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Fig. 10. — Rotation Measure Spread Function for a typical set of 
channels in the WSRT 350 MHz band. The RMSF is the normal- 
ized (unitless), complex response to polarized emission in Faraday 
space for a given set of A 2 sampling. Roughly 400 channels were 
used to construct this RMSF. 
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Fig. 11. — Rotation Measure Spread Function for a typical set of 
NVSS + 8 WSRT IF measurements. The RMSF is the normalized 
(unitless), complex response to polarized emission in Faraday space 
for a given set of A 2 sampling. 

iFrick et al.l (|2010T ). the method of RM Synthesis/ Clean 
has difficulty reproducing the correct phase information 
in the presence of multiple RM components. For this 
reason, we have neglected x m the solutions by fitting to 
p only. To exclude possible residual instrumental RM, 
which manifest near ±42 rad/m 2 due to the 17 M Hz 
modulation investigated bv lPopping fc Braunl (|2008f ). we 
searched for components in the range |RM| < 40 rad/m 2 . 
The RM Synthesis/Clean solutions for each source are 
summarized in Tables l6l - [T2l 

C. 2 component models for q(\ 2 ), u(X 2 ). This fit in- 
volves six parameters, with the amplitude of the frac- 
tional polarization, po, intrinsic polarization angle, xo ; 
and RM to be determined for each of two components. 
Because we expected (and sometimes found) multiple 
minima in x 2 m this six dimensional space, we mini- 
mized x 2 through a direct search of parameter space. 
The explored ranges were tailored somewhat to the indi- 
vidual sources, but typical values were polarized fraction 
(0, 0.1), RM (-25, 25) rad/m 2 , and xo (0, 180) degrees. 
The values presented in Tables [HI - 1 1 21 represent the min- 
imum of x 2 over these ranges. Note that there are no 
n7r ambiguities when fits are done in q,u space. A slice 
through this % 2 surface for the two RM parameters for 
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Fig. 12. — A cut through the \ 2 surface for our q, u vs. A 2 model 
grid search for 3C33S. The contours show the deepest minima in 
the surface, with the best fit RMs near -3 rad/m 2 (stronger com- 
ponent) and rad/m 2 (weaker component). Contour levels are 
at the probabilities of 10~ 6 , 10 -5 , 10 -4 , 10 -3 , 10 -2 , 2.5 X 1CT 2 , 
5 X 10~ 6 , HT 1 . 

3C33S is shown as an example in Figure Q21 Each value 
in this space represents the minimum value of \ 2 for fixed 
values of the two RMs, with all other parameters allowed 
to float. 

The errors in RM were calculated by normalizing Xmin' 
defining Xmin = dof. We then found the range of each 
of the two RMs for which the value of the normalized 
X 2 < (dof + 1), allowing the other five parameters to 
float. In a number of cases, there were additional minima 
within the x 2 < (dof + 1) range, so no errors are quoted 
and these RM values are shown in brackets. This proce- 
dure, of determining errors by adding 1 to the ySL- has 
a long history in the astrophysical literature (e.g., Avnil 
[19^ TWalil [19961) . but has very serious problems as dis- 
cussed below. The probability contour levels in Figure 
[12] were assi gned using the y 2 dis tribution for one degree 
of freedom (jLampton et al.lll976h as is appropriate when 
assigning errors to each individual RM, and allowing the 
other RM and all other parameters to float. 

D. Foreground rotation and depolarizing screen models 
forq(X 2 ), u(\ 2 ). We followed a procedure similar to that 
of the two component model, finding the minimum x 2 for 
the three parameter function 

p(A 2 )=<z(A 2 ^(A 2 ) 

= Po cxp(-24 M A 4 ) exp[2i( Xo + A 2 RM)] (14) 

similar to that described by I Burn! (Qj966). Errors in RM 
were determined in the same way as the two component 
model. 

We quote errors on RM using methods similar to those 
in the literature, so that our uncertainties can be com- 
pared to them. However, it is very rare in the literature 
to find x 2 values quoted for the fits, and therefore diffi- 
cult to evaluate whether the models used are appropriate 
or not. As we will show below, it is possible to get quite 



robust x(^ 2 ) °c A 2 behavior given an apparent RM quite 
different from the actual RMs for two component mod- 
els. Therefore, RM determinations using only x(A 2 ) and 
ignoring the fractional polarization behavior can provide 
no guidance regarding the appropriateness of the fit. 

The xl values in Tables El - HI are often much greater 
than unity, showing that these models are not an ad- 
equate representation of the data. In that case, the 
meaning of our errors is unclear. Our normalization 
°f Xmin = dof produces much more conservative er- 
rors th an simply adding 1 to Xmin- However, as pointed 
out bv lLampton et all () 19760 . this method produces the 
equivalent of a ratio of variances distribution, which 
has a very different probability distribution than x 2 it- 
self. In particular, they say "We stress again that if 
(Xmin dof), no formalism which uses distributions 
describing random fluctuations can provide the proper er- 
ror estimator. " Given this, our errors must be accepted 
only in the sense of providing comparisons to the lit- 
erature, and our recommendations for future work are 
described in Section |4"31 

3.2.4. Model Results and Comparisons 

The results of the various RM determinations for each 
source are shown in Tables [6] - [12] and in Figures [13] - HH 
Oscillation visible in the restored p(4>) profile (e.g., Fig- 
ure [TS]) arises from the sinusoidal residuals in q(4>) and 
u(4>), which are added to the clean components, and pro- 
vides a measure of the noise level in p((f>). In our tests, 
more aggressive cleaning reduced the oscillation in p by 
placing power in the clean components randomly across 
(p (thus reducing the residuals and producing nicer plots), 
but did not significantly change the amplitudes or Fara- 
day depths of the major RM components as reported 
for each source. If we look at the RM of the strongest 
component, we find that different models yield consis- 
tent results for some but not all sources. We give a 
brief discussion of the modeling results for each source 
here, and a more extensive discussion of 3C33S in the 
following subsection. One goal of this investigation is to 
explore the effect that ignoring the effects of depolariza- 
tion, including the presence of multiple strong RM com- 
ponents, may have on the findings of the traditional lin- 
ea r y(A ) meth o d. W e compare our findings with those 
of iTavlor et all £2009) for 

sources where RMs were re- 
ported in their study. Where appropriate, we use data 
from the VLA FIRS T (Faint Images of t he Radio Sky at 
Twenty centimeters; iBecker et al.l 119951 ) survey to sup- 
plement our analysis. 

NVSS J010616+125116. This source is resolved 
as a double source (separation ^60") in the original 
NVSS image, but appears as a single source when con- 
volved to the WSRT field resolution (325"x70"). We 
adopt the name of the brighter NVSS source (peak 
Ii. 4=102 mJy/beam); the secondary source is NVSS 
J010615+124210 (peak Ii. 4 =72 mJy/beam); the two 
sources have similar P1.4. The dominant RM is found 
near -9 rad/m 2 for the linear x(A 2 ), depolarizing screen, 
and RM Synth/Clean methods. The two component 
model, however, finds the dominant RM component near 
-5 rad/m 2 . It is possible that the relatively strong sec- 
ondary RM component found near +2 rad/m 2 in the 
two component fit has drawn the other solutions away 
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TABLE 6 

Modeling results for NVSS J010616+125116 



Model 


RMi (err) 


P01 


X01 


RM 2 (err) 


P02 


X02 


2 

"RM 


xi 




(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 




Linear x(A 2 ) 


-9.6 (0.12) 




80 










1.2 


Screen 


-9.5 (5) 


10.5 


78 








1.2 


10.2 


RM Synth/Clean 


-9.1 (0.04) 


3.1 




11.5 (0.2) 


0.5 








Two Component 1 


[-5.0] 


4.0 


75 


[2.0] 


2.2 


65 




8.8 



Brackets indicate multiple minima in x 2 surface - no RM error reported. 



TABLE 7 

Modeling results for 3C33S 



Model 


RMi (err) 


P01 


Xoi 


RM 2 (err) 


P02 


X02 


2 

a RM 


xi 




(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 




Linear x(^ 2 ) 


-6.8 (0.17) 




80 










3.9 


Screen 


-7.0 (0.15) 


8.5 


86 








1.0 


6.0 


RM Synth/Clean 


-6.7 (0.06) 


2.6 




8.0 (0.1) 


0.9 








Two Component 


-2.9 (0.1) 


6.7 


85 


-0.05 (0.2) 


4.8 


49 




2.1 



TABLE 8 

Modeling results for NVSS J011136+132437 



Model 


RMi (err) 
(rad/m ) 


Poi 
(%) 


Xoi 

(°) 


RM 2 (err) 
(rad/m 2 ) 


P02 
(%) 


X02 
(°) 


2 

(rad/m ) 


xi 


Linear x(-^ 2 ) 


-10.8 (0.31) 




130 










6.5 


Screen 


-10.75 (2.5) 


5.5 


130 








1.2 


9.3 


RM Synth/Clean 


-13.1 (0.1) 


2.0 




26.0 (0.3) 


0.6 








Two Component 


-11.2 (0.5) 


1.5 


140 


-24.2 (0.6) 


1.5 


175 




7.9 



TABLE 9 

Modeling results for NVSS J011204+124118 



Model 


RMi (err) 


P01 


Xoi 


RM 2 (err) 


P02 


X02 


2 

a RM 


xi 




(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 




Linear x(-^ 2 ) 


20 (1.9) 




60 










12.0 


Screen 


19.5 (0.5) 


7.5 


94 








1.0 


1.9 


RM Synth/Clean 1 


-32.8 (0.3) 


2.1 




34.4 (0.5) 


1.7 








Two Component 2 


[19.5] 


2.7 


95 


[-2.0] 


1.5 


160 




2.12 



1 Four strong features exist in the cleaned FDF, including one near 17 rad/m 2 . 
2 Brackets indicate multiple minima in x 2 surface - no RM error reported. 



TABLE 10 

Modeling results for NVSS J125630+270108 



Model 


RMi (err) 


Poi 


Xoi 


RM 2 (err) 


P02 


X02 


°RM 


xi 




(rad/m ) 


(%) 


(°) 


(rad/m ) 


(%) 


(°) 


(rad/m ) 




Linear x(-^ 2 ) 


4.6 (0.31) 




-175 










0.9 


Screen 


4.5 (0.25) 


1.0 


4 








0.0 


1.03 


RM Synth/Clean 


4.8 (0.3) 


1.1 




-17.7 (1.3) 


0.2 








Two Component 


4.5 (0.25) 


1.5 


10 


-4.0 


0.5 


130 




0.851 



TABLE 11 

Modeling results for NVSS J162408+605400 



Model 


RMi (err) 


Poi 


Xoi 


RM 2 (err) 


P02 


X02 


2 

a RM 


xi 




(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 




Linear x(-^ 2 ) 


-16.8 (0.4) 




65 










3.2 


Screen 


-17 (0.15) 


5.5 


64 








1.4 


1.2 


RM Synth/Clean 


-14.6 (1.5) 


0.5 




10.3 (0.4) 


0.2 








Two Component 1 


[-17.0] 


2.5 


95 


[-18.0] 


2.5 


55 




1.99 



Brackets indicate multiple minima in x 2 surface - no RM error reported. 
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TABLE 12 

Modeling results for NVSS J162740+603900 



Model 


RMi (err) 


P01 


X01 


RM 2 (err) 


P02 


X02 


4 

RM 


xi 




(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 


(%) 


(°) 


(rad/m 2 ) 




Linear x(-^ 2 ) 


-7.8 (1.5) 




65 










4.0 


Scrccn 1 


[-7.0] 


2.0 


34 








1.0 


1.8 


RM Synth/Clean 2 


-6.4 (0.8) 


0.5 




14.6 (0.5) 


0.4 








Two Component 


4.3 (0.5) 


1.0 


133 


15.0 (0.5) 


1.0 


160 




1.58 



1 Brackets indicate multiple minima in \ 2 surface — no RM error reported. 
2 Three strong RM components are present in FDF, including one near 4 rad/m 2 . 
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from the true intrinsic Faraday structure. The presence 
of multiple minima in the x 2 surface, however, casts un- 
cert ainty on the t wo co mponent result. For compari- 
son, iTavlor et all (|2009fl determined the RM of NVSS 
J010616+125116to be -16.8 ± 14.7 rad/m 2 ; no RM was 
reported for NVSS J010615+124210. 

3C33S. This source is also known as NVSS 
J010850+131831. The dominant RM found by the linear 
x(A 2 ), depolarizing screen, and RM Synth/Clean meth- 
ods are near -7 rad/m 2 . This is in disagreement with the 
two component modeling, which finds no significant com- 
ponent near -7 rad/m 2 ; rather, the dominant component 
is found near -3 rad/m 2 with a relatively strong second 
component near rad/m 2 . For comparison. ITavlor et all 
(|2009ft d etermined the RM to be 3.4 ± 1.9 rad/m 2 . In 
addition, lLaw et aT] ([201 ll ) performed RM Synthesis on 
3C33S using two bands, each 100 MHz wide, centered at 
1.43 and 2.01 GHz with the ATA. After cleaning they 
found a single RM at -12.3 ± 0.4 rad/m 2 . That they 
found a single RM is not unexpected, considering the 
FWHM of their RMSF of 141 rad/m 2 , but the RM value 
found would not fit our 350 MHz observations. Given the 
high signal to noise in our Q, U data, this is the strongest 
case yet for interference between two strong RM com- 
ponents causing other methods to misinterpret the true 
Faraday structure. In the next section, we will use ide- 
alized models to demonstrate how this comes about. 

NVSS J011136+132437. The dominant RM compo- 
nent is found near -11 rad/m 2 for the linear x(A 2 ) an d 
depolarizing screen models. The two component method 
finds equal amplitudes for both RM components, with 
one near -11 rad/m 2 and the other near -24 rad/m 2 . 
It doesn't appear that a secondary component has af- 
fected the outcome of the single component methods. 
RM Synth/ Clean nearly agrees, finding the dominant 
RM component at -13 rad/m 2 . The secondary RM com- 
ponent found by RM Synth/Clean and the two compo- 
nent model are in disagreement, however, in both loca- 
tion and relative amplitude. All methods have a high 
xl, suggesti ng that no solution is to be trusted. For 
comparison, ITavlor et all ([20091 ) determined the RM to 
be — 13. 8 ± 3.3 rad/m 2 , in agreement with our findings. 

NVSS J011204+124118. The dominant RM is found 
near +20 rad/m 2 for the linear x(A 2 ), depolarizing 
screen, and two component models. The two compo- 
nent fit finds a secondary component with P02/P01 > 0.5 
near -2 rad/m 2 , but it doesn't appear to have affected 
the outcome of the single RM methods. RM Clean finds 
the dominant RM component near -33 rad/m 2 , but three 
other peaks of significant amplitude are found in the 
Faraday spectrum, including relati yely strong c ompo - 
nents near +34 and +17 rad/m 2 . ITavlor et all (|2009l ) 
do not report a RM for this source. 

NVSS J '125630+270108. All four methods find the 
dominant RM component to lie near +4.5 rad/m 2 . Sec- 
ondary components for RM Synthesis and the two com- 
ponent method are of relatively weak amplitude, and 
likely do not contribute significantly to the solutions 
found by the single RM methods. Due to the lack of 
depolarization from 1.4 GHz to 350 MHz, it is not sur- 
prising that the traditional linea r fit to \(X 2 ) is in ag ree- 
ment with the other methods. ITavlor et al.l (|2009f) do 
not report a RM for this particular source (unresolved in 



both NVSS and FIRST), but using their data we deter- 
mined the weighted mean RM of the 1 7 sources within 
2° (with an entry in ITavlor etall 120091 ) to be «2.5 ± 1 
rad/m 2 . This suggests that a Galactic foreground (rotat- 
ing) screen is the single dominant component of Faraday 
structure for this source, a situation which is unique in 
our modeling results. 

NVSS J162408+ 605400. The dominant RM compo- 
nent is found to lie near -17 rad/m 2 for the linear x(A 2 ) 
and depolarizing screen models, while RM Synth/Clean 
finds the dominant component near -15 rad/m 2 . The two 
component model shows two equal amplitude RM com- 
ponents near -17 rad/m 2 and -18 rad/m 2 . The presence 
of multiple minima in the x 2 s urface casts uncerta inty 
on the two component solution. ITavlor et al.l (|2009f ) do 
not report a RM for this source. The oscillation visible 
in the cleaned FDF shown in Figure [18] is due to the 
low amplitude of the clean components relative to the 
amplitude of the residuals in RM Clean. As mentioned 
in the beginning of this section, cleaning further would 
reduce the level of apparent oscillation in the cleaned 
FDF, but would not appreciably change the locations or 
amplitudes of the fitted RM components. 

NVSS J162740+603900. The dominant RM compo- 
nent is found by the linear x(A 2 ), depolarizing screen, 
and RM Synth/Clean methods to be near -7 rad/m 2 . 
The cleaned FDF displays three strong RM features, near 
-7, +4, and +15 rad/m 2 , but it is likely that the compo- 
nents at -7 and +4 rad/m 2 are blended, contributing 
power to each other and increasing their peak ampli- 
tudes. These results contrast with the two component 
method, which finds two dominant RM components of 
equal amplitude near +4.5 and +15 rad/m 2 . Again, it 
seems likely that two RM components are interfering in 
a way which confounds the other methods. ITavlor et al.l 
(2009J) do not report a RM for this source. 

A comparison between the different methods of deter- 
mining RMs for each source are shown in Figure [20] As 
expected, the linear x(A 2 ) an d depolarizing screen fits 
agree well for the dominant RM value since nir angle 
shifts were inserted into the data for the x(A 2 ) fits to 
best match the depolarizing screen models. The RM 
of the dominant component found by RM Synth/Clean 
agrees fairly well with the linear x(A 2 ) fit method for six 
of the seven sources, although only three agree within 
the formal errors. The dominant RM found by the two 
component model fit, however, finds agreement with the 
linear x(A 2 ) fit method in only four of the seven sources 
analyzed. It is apparent that in six of the seven sources 
the traditional linear fit to x(A 2 ) is incapable of provid- 
ing a description of the source's true Faraday structure, 
instead providing what may be referred to as a "charac- 
teristic" RM. This is due to the fact that fitting to x(A 2 ) 
does not consider the behavior of p(A 2 ), which is variable 
in most of our sources when the measurements across a 
large range of A 2 are considered. One must consider 
depolarization models, such as the depolarizing screen 
or interference between multiple RM components, if the 
true Faraday structure is to be described. 

We note that the x 2 values for these fits are quite 
high in many cases, suggesting more complicated mod- 
els would be needed to properly fit the data. Some of 
the data appear anomalous when the apparent behavior 
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of neighboring points is taken into account. These data 
could be contaminated by residual instrumental prob- 
lems; we have attempted to incorporate these effects into 
our errors. Given infinite resources, our instrumental er- 
rors would approach zero; polarization calibration is no- 
toriously difficult at low frequencies and we must there- 
fore proceed with our best effort, given the current tech- 
nological limitations. By removing "anomalous" data 
we would be biasing the solutions toward simple Fara- 
day structures in the model fitting and degrading the 
ability of RM Synthesis to resolve multiple components 
closely spaced in Faraday depth. In addition, we found 
that some of the discrepancies between the fits were due 
to flaws in the techniques themselves, which we discuss 
in Section 1331 using the case of 3C33S. 

3.3. 3C33 South 

3C33S is the southern lobe of the z = 0.059 radio 
galaxy 3C33 near the pointing center of the Aries-Pisces 
field. At an off-axis radius of ~650", its mean p of ~4% 
in the 350 MHz band is well above the mean WSRT in- 
strumental contribution of <0.5% at that radius. In the 
NVSS image, convolved to the common beamsize used 
for our WSRT images, 3C33S displays a fractional po- 
lar ization of 10% . Prio r studies at A6 cm and A2 cm 
bv lRudnick et al.l (|198lD . which resolve the structure of 
the lobe, find the fractional polarization ranging from 
~15% at the radio peak to more than ^60% in the 
lower surface brightness regions. Previous studies have 
quoted an integrated RM of -12 rad/m 2 for 3C33 (e.g . , 
iBerge fc Seielstadl[T967llSimard-Normandin et al.lll98ll ). 
although this RM det ermination may be contaminated 
by the northern lobe. iRudnickl (|1988l) finds the RM to 
be ss-7 rad/m 2 in the southern hotspot region using un- 
published 20 cm and 6 cm data. 

Table [7] and Figure [T4l summarize the modeling results 
for 3C33S. The model fit for the depolarizing screen also 
yielded RM = -7 ± 0.15 rad/m 2 , as did the x(A 2 ) fit 
alone (RM = -6.8±0.17 rad/m 2 ). Similarly, RM Synthe- 
sis/Clean found a dominant component with p = 2.5% at 
-6.7 rad/m 2 , with a weaker p = 0.6% component at +8 
rad/m 2 . This result is robust for various weighting of the 
NVSS sample with respect to the WSRT samples before 
computing the FDF. While increasing the weighting of 
the NVSS sample can have a large effect on the sidelobe 
level and structure in the RMSF and, hence, the con- 
structed FDF, RM Clean yields a similar solution each 
time: a dominant peak near -7 rad/m 2 and a secondary 
peak near +8 rad/m 2 . Thus, there appears to be good 
agreement between the literature, depolarizing screen, 
x(A 2 ) and RM Synthesis/Clean results that the domi- 
nant RM component in 3C33S is at -7 rad/m 2 . Weaker 
RM components in the FDF, such as the p = 0.4% one 
at +23 rad/m 2 , are increasingly unreliable (see Figure 

QU). 

However, our two component fit to these same data 
give quite different values, -3 and rad/m 2 . Which of 
these determinations is correct? While direct comparison 
of the xt values for the best model fits is inappropriate, 
it is clear from Figure [TJ] that the two component model 
provides a better fit to the data than the depolarizing 
screen model. It does a much better job of explaining 
the two longest wavelength observations (particularly in 
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Fig. 13.— Polarization diagnostics for NVSS J010616+125116. 
Model fits (lines) are plotted over the observed data (points). Top 
panel: two component model. Middle panel: depolarizing screen. 
Bottom panel: magnitude of the cleaned fractional FDF (solid line) 
and rms of the residuals (horizontal dashed line). The linear x(-^ 2 ) 
fit is omitted since it is nearly identical to the depolarizing screen. 

q), where p rises from a minimum near A 2 ~ 0.8 m 2 . The 
observation of a minimum in p is compelling evidence 
against a simple depolarizing screen. In any case, a slight 
change in the errors assigned to the original data points 
could change the relative goodness of fit f or these two al- 
terna tives. Polarization data at A = 9 cm (Rudnick et al. 

agrees in p with our two component and depolariz- 
ing screen models, but have been excluded from the fits 
since they were integrated over both lobes and therefore 
not reliable for these purposes. 

We note that such a discrepancy can be quite impor- 
tant depending on the scientific issues under investiga- 
tion. First, these two models (a single component at -7 
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Fig. 14. — Polarization diagnostics for 3C33S. Same layout as 
Figure [13] 



0.04 
0.02 


-0.02 
-0.04 
-0.06 
-0.08 




0.2 0.4 0.6 0.8 

% 2 (m 2 ) 



2 Component Model 

0.04 : 
\l \ 0.02 - 


' -0.02 
-0.04 
-0.06 
-0.08 



0.2 0.4 0.6 0.8 

% 2 (m 2 ) 








♦ 





-0.06 0.06 



Depolarizing Screen Model 




0.06 0.06 




-20 20 40 

Faraday Depth (rad/m 2 ) 



Fig. 15.— Polarization diagnostics for NVSS J011136+132437. 
Same layout as Figure [13] 



rad/m 2 or two components, at -3 and rad/m 2 ) repre- 
sent quite different physicai structures in the source. For 
exampie, the magnetic field in 3C33 S very closely t racks 
the bow-shocked shape leading edge (Rudnick 1988;), and 
a small toroidal sheath could give rise to two domi- 
nant RM components. Alternatively, the surrounding 
medium might have a depolarizing screen with very fine 
scale structure (^Cl", ~1 kpc) that is independent of 
the geometry of the source. If we were not interested 
in the Faraday structure, we could simply look at the 
weighted mean of the two component fit, which yields 
-1.7 rad/m 2 . However, the difference between this value 
and the -7 rad/m 2 from other models represents a fac- 
tor of greater than 4 in any derived densities or mag- 
netic field strengths. If similar discrepancies are found 
at shorter wavelengths, e.g., 1 GHz, then the physical 



parameters involved would be ~10x larger. 

The discrepancy between the -7 rad/m 2 and weighted 
-1.7 rad/m 2 fits is not due to the inaccuracies of the 
measurements, as determined from the formal errors. 
The error in our x(A 2 ) fit is small (RM = -6.8 ± 0.17) 
and the errors in the two component fit are even smaller 
(RM = -2.9 ± 0.1). Thus, using our linear x(A 2 ) fit we 
would have ruled out the two component weighted mean 
with high confidence. Similarly, our RM Synth/Clean 
results would have ruled out the two component fit. It 
might be further argued that we shouldn't have expected 
to distinguish between values of -7 and -1.7 using this 
method, since the FWHM of the RMSF's main lobe is 
~12 rad/m 2 . This argument ignores the standard prac- 
tice of quoting uncertainties in the location of a peak 
at a value of ~FWHM/(2xsignal:noise). In the case of 
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Fig. 16.— Polarization diagnostics for NVSS J011204+124118. 
Same layout as Figure 1131 

3C33S, the error in the dominant RM peak (from Gaus- 
sian fitting to the cleaned FDF) is 0.06 rad/m 2 . Again, 
we would have ruled out the (unresolved) combination of 
peaks near -3, rad/m 2 with high confidence. 

Because the discrepancies between the results of var- 
ious models, in particular RM Synthesis/Clean, were 
much larger than our calculated errors, we carried out 
a series of experiments with infinite signal to noise mod- 
els using two RM components. 

4. EXPERIMENTS WITH TWO RM COMPONENTS 

4.1. RM Synthesis 

Although many different polarization diagnostic ex- 
periments could (and should) be done, we focused on 
two-component models for several reasons. First, a two- 
component model produced a good fit to the 3C33S data. 
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Fig. 17.— Polarization diagnostics for NVSS J125630+270108. 
Same layout as Figure 1131 



Second, two Faraday components might be a reasonable 
expectation for double radio galaxies that are unresolved. 
In addition, when angular resolution becomes sufficient 
to resolve depolarizing Faraday screens, there will always 
be places where the beam overlaps two neighboring struc- 
tures, produ cing two Faraday components. Finally, re- 
cent work by lLaw eTaTl (|20Tll) . where RM Synth/Clean 
was performed on 37 polarized radio sources using the 
Allen Telescope Array (ATA), showed that two or more 
components were detected with high confidence in ~25% 
of their sources. We now discuss a few simple experi- 
ments to demonstrate some of the potential pitfalls when 
RM Synthesis/Clean is employed. 

Our first experiment was to adopt a model fixed to the 
best two-component fit to the 3C33S data, with compo- 
nents at -3 and rad/m 2 . Synthetic q,u spectra were 
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Fig. 18.— Polarization diagnostics for NVSS J162408+605400. 
Same layout as Figure [13] 
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Fig. 19.— Polarization diagnostics for NVSS J162740+603900. 
Same layout as Figure [13] 



constructed for the same NVSS + 8-WSRT A 2 loca- 
tions as in our previously discussed observations. The 
results of RM Synth/Clean are shown in Figure [2T1 It 
bears a remarkable resemblance to the observed FDF for 
3C33S, displaying a single dominant peak near RM ~ —7 
rad / m 2 and a low amplitude secondary feature near RM 
~ +8 rad/m 2 , even though the input RMs were at -3 and 
rad/m 2 . Thus, the FDF, with or without cleaning, pro- 
duces RM power at what we can now state is the wrong 
value, since we know the input model parameters. This 
is true whether you examine the clean components at 
high RM resolution or their convolved version which re- 
flects more closely the limitations in resolving multiple 
RM components. In the convolved case, one would ex- 
pect the FDF to still reflect the weighted mean of the 
input components; it does not. 



Another case, demonstrating the impact of the relative 
phase of the two polarized components, involves using 
two components of equal amplitude with RMs of -15 and 
rad/m 2 . These are separated by more than the FWHM 
of the RMSF, 12 rad/m 2 , constructed from ^400 chan- 
nels in the WSRT 350 MHz band and shown in Figure 
[lT)l Nominally, then, they should appear well-separated 
in the FDF. Figure |2"21 shows the results of using four dif- 
ferent values for the difference in \q for the two compo- 
nents. In three cases, RM Synthesis/Clean successfully 
resolved the two components. In the fourth case, with a 
difference in xo of 90°, the raw FDF was dominated by 
a single peak near the mean RM of -7.5 rad/m 2 , along 
with considerable sidelobe power. Cleaning produced an 
apparent triple component structure, with power at RMs 
of -17, +2 and -7.5 rad/m 2 , instead of the input values 
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Fig. 20. — Comparison of the methods for RM determination for the seven modeled sources. Upper left: depolarizing screen vs. x(-^ 2 ) 
fit. Upper right: RM Synth/Clean vs. x(\ 2 ) fit. Lower left: two component fit vs. x(^ 2 ) fit- The two strongest RMs are plotted for RM 
Synthesis/Clean and two-component model fitting, connected by thick dashed lines for each source. Pointsize for the primary component 
is fixed, while the pointsize of the secondary co mpon ent (relative to that of the primary) is proportional to the ratio of amplitudes for the 

6112 Lower right: RM Synth/Clean vs. two component fit. Large and small points show 
Errors from the fitting techniques are plotted, but are smaller than the pointsize for 



RM components (i.e. P2/P1), as listed in Tables 
RMi and RM2, respectively, as listed in Tables 



6TT2 



most sources. 



of -15 and rad/m 2 . 

We have also asked colleagues to run these and other 
models through their own RM Synthesis and cleaning 
programs, to verify that coding problems were not at 
fault. The above problems with RM Synthesis and Clean 
are robust to its exact implementation. They occur when 
the separation of the RMs is on the order of the FWHM 
of the RMSF, 2y / 3/AA 2 . Under these conditions, the 
number of cycles of Q and U within the bandwidth differ 
by one or less for the two components. RM synthesis is 
therefore not able to reliably resolve them into separate 
Fourier components. However, the two components do 
not simply blend in this case, as two nearby sources would 



blend in total intensity. Instead, they interfere to create 
complicated structures in Q(X 2 ), U (X 2 ) (i.e. P(\ 2 ) and 
x(X )) which causes RM Synthesis to put power at values 
other than the input RM. In some cases, RM Clean is able 
to recover from this interference; in other cases it is not. 

Some of the shortcomings of RM Synthesis arise not 
from a fault in the technique, but rather a limit of our 
measurement abilities. One can use radio aperture syn- 
thesis as an analogy from which to draw insight; limita- 
tions in baseline sampling for aperture synthesis are in 
some ways analogous to limitations in A 2 sampling for 
RM Synthesis. However, the RM interference that we 
have illustrated here is considerably more complicated. 
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Fig. 21. — Cleaned FDF for the best fit two-component model 
of 3C33S. The two input RM components are at -2.9 and -0.05 
rad/m 2 as discussed in the text, but the dominant peak in the 
Faraday spectrum is near -7 rad/m 2 with a secondary feature 
near +8 rad/m 2 . Vertical dashed lines show the location of the 
two input RMs, -2.9 and -0.05 rad/m 2 . 

These experiments reflect the interference between two 
RM components and are reminiscent of other types of 
interfer ence that are better understo od. Polarization 
canals (jShukurov &; Berkhuijsenll2003t ). e.g., do not rep- 
resent actual dips in polarization, but simply the interfer- 
ence, in one beam, between two components separated by 
90° in polarization angle at some observed wavelength. 
Similarly, rotation measure involves the trend of x(A 2 ) 
over a range of wavelengths, and the mapping between 
multiple RM components and x(A 2 ) is not yet fully un- 
derstood. This illustrates the need for sufficiently broad 
A 2 coverage in polarization observations when perform- 
ing RM experiments, where detection of potential max- 
ima or minima in p(A 2 ) is also critical to help diagnose 
the Faraday structure. These methods are also subject 
to a variety of degeneracies, some of which we illustrate 
in the following section. 

4.2. Pseudo-X 2 behavior 

Another insidious quality of two component models is 
that they commonly produce 

RM{\ 2 ) = dx{\ 2 )/d\ 2 w constant (15) 

over substantial ranges in A 2 space. Although it may 
be obvious that sparsely sampled data (especially using 
only two or three A 2 data points) could lead to mistakes, 
it is assumed that continuous sampling over a significant 
range of wavelengths (e.g., (A 2 raa2 ,-A 2 nm )/A 2 > 0.25) can 
verify whether RM(A 2 ) ps constant. This is not always 
true, as we now illustrate. 

Figure [53] shows five different models, all of which 
produce excellent RM(A 2 ) s» constant over the WSRT 
350 MHz band, which covers ~35% in A 2 space. In addi- 
tion, three of these models would also yield the same ex- 
cellent A 2 behavior with an additional point at 1.4 GHz 
(e.g., NVSS). The model parameters are listed in Table 
1131 If one were examining the behavior of x(A 2 ) alone, 
as is done in most of the existing literature, there are a 
wide variety of two component models which easily fit 
the data but have very different values of RM than the 
one observed. The key to ruling out such two component 
models, and thus to have a reliable determination of RM, 



TABLE 13 

Model Parameters for Long Wavelength 
Pseudo-A 2 Experiment 1 



Model 


RMi 


RM 2 


Pi 


P2 


XI 


X2 


ID 


(rad/m 2 ) 


(rad/m 2 ) 


(%) 


(%) 


(°) 


(°) 


Modi 


1 




1 





90 




Mod2 


-0.5 


2 


0.66 


1 


-43 


-123 


Mod3 


-4.4 


6.3 


1 


1 


-30 


-145 


Mod4 


1 


2 


1.1 


.55 


-75 


3 


Mod5 


-97 


99 


1 


1 


-35 


35 



1 See Figure E3 



TABLE 14 

Model Parameters for Short Wavelength 
Pseudo-A 2 Experiment 1 



Model 


RMi 


RM 2 


Pi 


P2 


XI 


X2 


ID 


(rad/m 2 ) 


(rad/m 2 ) 


(%) 


(%) 


(°) 


(°) 


SI 


1110 


750 


0.7 


0.35 


-92 


92 


S2 


1650 


500 


0.5 


0.55 


-110 


102 


S3 


1400 


-50 


0.665 


0.35 


76 


125 


S4 


1000 




1 





88.85 





1 See Figure [24] 



lies in their p(A 2 ) behavior, which is quite different for 
each model. A better way to avoid these mistakes is to 
simply fit the function p = pe 2lx = q + iu to the q{\ 2 ) 
and m(A 2 ) data, and determine whether a satisfactory fit 
has been achieved. 

It is tempting to assume that there is a "short wave- 
length" limit where these problems can be safely ignored. 
We now show that is not true. First, we define a "short 
wavelength" set of observations as one in which there 
is reasonable sampling in A 2 space and Aj^-C Aj, ra . 
Thus, one can verify whether RM(A 2 ) « constant down 
to effectively zero wavelength. The models described in 
Table [T4l and shown in Figure |2~41 demonstrate that this 
does not exclude two component models with RMs very 
different than the ones measured by fitting RM(A 2 ) = 
constant. In the examples shown, the x(A 2 ) data alone 
follow very closely a constant RM = 1000 rad/m 2 . How- 
ever, they actually contain components that range from 
-50 to 1650 rad/m 2 . Again, the key is to examine the 
p(A 2 ) behavior, as seen in Figure IM1 or better, as noted 
before, to fit a model directly to the q(X 2 ) and u(A 2 ) 
data. 

We also show in Figure [Ml the full RM(A 2 ) for these 
models. The wide variations in this number show that 
when there are two interfering components, measuring 
the RM using data at closely spaced wavelengths, or only 
using very sparse sampling, can render the observed RM 
virtually meaningless. 

4.3. Recommendations regarding RM measurements 

There is no simple prescription for producing reliable 
rotation measures because it depends on the specific sci- 
entific goals. We begin the discussion of those issues 
below, but here we simply offer some general guidelines 
to inform the future practice of Faraday structure deter- 
minations: 

• Fitting of models to q(\ 2 ) and w(A 2 ) (or equiva- 
lently, p(A 2 ) and x(A 2 )) is the only reliable way 
to determine the underlying Faraday structure. In 
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Fig. 22. — Illustration of the effect of relative phase between two RM components upon the results of RM Synthesis/Clean for various Axo 
configurations. Although the two co mpo nents are separated by more than the FWHM (12 rad/m 2 ) of the RMSF, RM Synthesis/Clean, 
using the same channels as in Figure ITDl fails to properly reproduce the solution for certain relative phases. The model RMs are at -15, 
rad/m 2 , shown by vertical dashed lines in the FDFs. Plotted in each panel are: top left: Fractional polarization, q (dashed), u (dotted), 
p (solid); top center: polarization angle (radians); top right: q vs. u; bottom left: Dirty FDF; bottom center: RM Clean clean 
components; bottom right: Cleaned FDF. 



particular, results derived from x(A 2 ) alone or RM 
Synthesis alone are subject to large ambiguities. 

• RM Synthesis/Clean, as it is currently imple- 
mented, can serve as a first order indicator of the lo- 
cation of power in Faraday depth space, and guide 
more detailed modeling. 

• Plots of q vs. u provide another useful diagnostic 
of the appropriateness of any models. 

• Results for RM determinations should always spec- 
ify not only the formal errors, but also the x 2 or 
RMS residuals of the fits. This, along with docu- 
menting the coverage in A 2 space, will allow for an 
analysis of what ambiguities are permitted by the 
data. 

• The allowable space for ambiguities can be signif- 
icantly reduced by broadening the A 2 coverage, 
increasing the sampling, and ensuring that regions 
of A 2 space are observed where RM(A 2 ) ^ constant 
and p(A 2 ) 7^ constant. 

• Scientifically useful results are possible in the pres- 
ence of ambiguities if the underlying assumptions 
are both documented and valid, as discussed in Sec- 
tion [5] 



• Alternative methods of parameter determination, 
such as Maximum Likelihood, should be considered 
in the presence of low signal to noise. In this case, 
least squares fitting may yield a low % 2 statistic, 
but may not necess arily yield the appropr iate solu- 
tion. For example, IGuidetti~ al. ( 2008) use four 
frequency samples in a linear fit to x(A 2 ) to deter- 
mine RMs for a number of cluster sources. In the 
limit of infinite signal to noise, it doesn't matter 
how closely spaced the points are; with four points 
and two parameters, i.e. dof=2, x 2 ~ 1 would 
truly signify a good fit. However, we note that 
for each source two of their samples are at nearly 
the same frequency, and these measurements agree 
within errors. This essentially guarantees a value 
for xt of order unity, perhaps giving false confi- 
dence in the appropriateness of the model. Thus, 
if minimization of \l lS to be used, we caution that 
the effective degrees of freedom should first be care- 
fully considered. 

5. DISCUSSION 

5.1. Depolarization 

The median depolarization ratio between 350 MHz and 
1.4 GHz for our seven modeled sources is P350/P1.4 ~ 0.3. 
This is the same as the median of the upper limits for 
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Fig. 23. — Various models illustrating the importance of consid- 
ering polarization amplitude as well as angle in the long wavelength 
regime. A linear fit to x(-^ 2 ) yields the same "characteristic" RM 
in each case, but inspection of the amplitude behavior reveals the 
co mpli cated nature of the various Faraday structures listed in Ta- 
ble [13] Top: Polarization angle vs. A 2 . Bottom: Polarization 
amplitude vs. A 2 . 

our sample as a whole. If this depolarization is due to 
a random foreground screen (and not to the interference 
between two components), then this corresponds to a 
Burn law grm ~ 1 rad/m 2 , where the depolarization 
is exp(— 2(7 2 ?M A 4 ). It is likely that the overall sample 
is even more depolarized, since we observed the upper 
limits to drop as the polarization fraction at 1.4 GHz in- 
creased (see Figured]). This result has two implications, 
one for observations at low frequencies and one concern- 
ing the environment of radio galaxies. 

Assuming a Faraday dispersion as above, we can es- 
timate, e.g., the depolarization that would be observed 
by LOFAR0 which has a high frequency band covering 
120-240 MHz, and a low frequency band covering 30- 
80 MHz. If the Burn law were to remain roughly accu- 
rate for integrated polarizations, then the depolarization 
would peak at 0.008 at the high end of the high frequency 
band, and drop by many orders of magnitude at low fre- 
quencies, essentially making polarizat i ons u ndetectable. 
However, as pointed out bv iTribblel (|1991h . the falloff 
from a Gaussian depolarizing screen is likely to be con- 
siderably slower, dominated by the small patches around 
extrema in RM, where the RM gradient is near zero. If 
we start with a characteristic integrated polarization of 
~3% at 1.4 GHz and extrapolate with only a A 2 depen- 
dence from our depolarization results at 350 MHz, then 

2 http://www.lofar.org 
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Fig. 24. — Various models illustrating the importance of consider- 
ing polarization amplitude as well as angle in the short wavelength 
regime. A linear fit to x(-^ 2 ) yields the same "characteristic" RM 
in each case, but inspection of the amplitude behavior reveals the 
co mplic ated nature of the various Faraday structures listed in Ta- 
ble I14I Top: Polarization angle vs. A . Center: Polarization 
amplitude vs. A 2 . Bottom: RM(A 2 ). 

we would expect fractional polarizations of 0.1% - 0.5% 
in LOFAR's high band, and 0.01% - 0.05% in the low 
band. These are not likely to be detectable. It is not 
clear, at present, whether even well-resolved extragalac- 
tic sources will have small enough Faraday dispersions to 
be observed in polarization at these low frequencies. 

Using the m o re p hysical units introduced by 
iGarrington et al.l (|1991l ). our observed characteristic 
lower limit to the Faraday dispersion is ^1.5 cm~ 3 /iG 
pc. For the purposes of calculating some very rough 
estimates of what these limits mean for field strengths 
around radio galaxies, we assume that the depolariza- 
tion occurs in a foreground screen completely unaffected 
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by the radio galaxy. Assume that we need ~10 inde- 
pendent patches across a 100 kpc source in order to de- 
polarize it, and a fiducial electron density of n e =10 -3 
cm~ 3 . The resulting magnetic field is then B///G > 
0.1(n e /10" 3 )(r/10)(cm 3 kpc -1 ), where r is a typical 
scale size of magnetic field fluctuations and we have ig- 
nored the y/N averaging along each line of sight for this 
order of magnitude calculation. Fields of this strength 
are less than those found in clusters of galaxies, but 
greater than expe cted in the mor e filamentary WHIM 
outside of clusters (|Rvu et al.ll2008D , especially if one fac- 
tors in the much lower densities in those regions. 

Thus, radio galaxies appear to be associated with ther- 
mal, magnetized plasmas with much higher values of 
n e B than expected for filamentary regions, but similar to 
those found in clusters. This could result because of the 
bias for radio galaxie s to be found in high density regions 
(|de Zotti et al.l [2010h . Alternatively, effects very local 
to the parent galax y, such as emission-line regions (e.g., 
IPedeltv et al.l IT 989) could be responsible for the ubiqui- 
tous depolarization. This leaves very few radio galaxies 
available to probe cosmologica l filaments, exc ept, per- 
haps some Mpc-scale sources (|Sari D allill2009l) . As we 
seek to understand the causes of the ~1.5 cm -3 /xG pc 
limits, however, it will also be important to readdress the 
questions of internal depolariz ation, e.g., due to a mix- 
ing layer ( Bicknell et al.l 1 19901 ) between the radio source 
(with possibly much higher fields) and its low density 
environment. 

5.2. Science Implications of RM Ambiguity 

We have shown that there is considerable ambiguity 
(sometimes >100%) in the determinations of RMs us- 
ing the methods universally used in the literature, and 
even in the more recent RM Synthesis technique. We 
now briefly examine the implications this has for differ- 
ent types of scientific investigations. 

5.2.1. Galactic Foreground 

The use of polarized extragalactic sources to charac- 
terize the magnet ic field structure of our Galaxy has a 
long history (e.g.. iSimard-Normandin fc KronbergllT980L 
iBrown fe Tavlorll200lD plus a maior recent advancement 
(jTavlor et al.ll2009l ). Our investigations do not reveal any 
(signed) bias in RM determinations, therefore, we would 
expect that the average RM of a group of extragalactic 
sources in some area of the sky should be a fair mea- 
sure of the true value. However, the structure function 
of galactic fluctuations will have contributions from the 
RM ambiguities discussed here, as well as attempts to 
measure the intrinsic differences in RM between sources, 
especially on the smallest scales. 

5.2.2. Fluctuations through Galaxy Clusters 

The situation with respect to cluster measurements 
is much more complicated. The clusters are expected 
to have fields that are tangled on scales substan- 
tially smaller than the cluster. Therefore the mean 
RM of a distant extragalactic source seen through the 
cluster should be zero, but the scatter in such RMs 
should be larger than f or background so u rces not seen 
through clusters (e.g.. IKim et al.l I1991L iClarke et al.l 
Bonafed e et al.ll2010D . In this case, the quantity 



being measured is the RM scatter, which will be in- 
creased because of the ambiguities discussed in this pa- 
per. In the ideal world, this scatter should be no differ- 
ent for sources seen through clusters ( "the sample" ) than 
for sources not seen through clusters ("the control"), so 
again the measurements should be unbiased. 

However, unless the sample and control have exactly 
the same properties, both intrinsically and in terms of 
the observations leading to their RMs, it is impossible to 
know how the RM ambiguities would affect their compar- 
ison. For example, if RM determinations include a short 
wavelength point for some sources, as opposed to oth- 
ers, a different range of possible underlying RMs will be 
present for the two cases. Or, if the sources in the sample 
or control are statistically different physically (e.g., FRI 
vs. FRII sources), then the ambiguities can have differ- 
ent effects and contaminate the test. All of t hese prob- 
lems a r e present in the w ell-cit ed studies by IKim et all 
(fl990T) . IKim et al.1 (fl99l. and IClarke et al.l (I2001D. as 
discussed by iRudnick fc Blundelll (|2004ft . Similar con- 
tamination can be present if RMs from one experiment 
are compared to RMs from another, with different wave- 
length cove rage, different editing for non-A 2 behavior, 
etc. (e.g., IJohnston-Hollitt fc Ekersl I2004D . Since our 
modelling shows that RMs can be affected by factors 
of order unity, it is not possible to assess the reliabil- 
ity of these cluster background experiments. Two types 
of studies are required to address this issue. First, the 
prevalence of multiple R M component s with in observing 
beams must be estimated; IPizzo et al.l (|2011f) found mul- 
tiple components at ~15" resolution in all three radio 
galaxies near the center of Abell 2255. Second, statisti- 
cal predictors are needed to quantify the likely errors in 
RM for a given distribution of multiple components. 

5.2.3. Faraday structure of radio galaxies 

Increasingly detailed studies of the Faraday struc- 
ture of individ ual radio galaxies ar e now becoming 
available ('e.g.. lLaing fc Bridle! 119871 iZavala fc Tavlorl 
20021 lLaing et all I2006L lO'Sullivan fc Gabuzdal I200S 



Govoni et al.l 120101 ) . In many cases, the rotation mea- 
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sures are assumed to be entirely in the unperturbed fore- 
ground, and thus a fair measure of the magnetic field 
structure of the environment, usually a clust e r of g alax- 
ies. However, as suggested by I Bicknell et al.l (|1990l ) and 
Rudnick & Blundcll (2003), and now demonstrated con- 
vincingly by lLaing et al.l (|2008D , the radio source itself 
may change the observed RM structure. This issue aside, 
the question remains how the newly described RM am- 
biguities could affect these measurements. 

The studies of individual radio galaxies involve higher 
order characterizations of the RM distribution, such as 
the structure function, so they are much more sensitive 
to possible ambiguities. In addition to increasing the 
overall scatter in RMs, contributions from ambiguities 
are likely to change as a function of scale. If we assume 
an unperturbed foreground screen, then when the ob- 
serving beam is much smaller than the smallest angular 
scale of RM variations, a single component dominates 
and the RM dete rmination can be free of ambiguities. 
iFeain et al.l (|2009l ) took advantage of this situation in a 
Faraday structure study of the radio lobes of Centaurus 
A, using background sources. 

In the limit where the observing beam is much larger 
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than the characteristic scale of variations, then we ap- 
proach the Burn limit of a depolarizing screen, and the 
effect o f ambigu i ties is minimized. However, as pointed 
out by iTribbld (|1991h . the situation is typically much 
more complicated, and the observed polarized emission 
will be dominated by regions where the angular RM dis- 
tribution is at an extremum, with only small gradients. 
The emission is then a complicated function of the beam 
size and the angular structure of the magnetic field fluc- 
tuations. Detailed modeling is required in such cases, 
and it is not certain whether a clear diagnosis of the 
Faraday structure is possible, in practice. In particular, 
we may not be able to distinguish between the physi- 
cally distinct cases of fully external screens, thin mixing 
layers of relat ivistic and thermal e mission or fully mixed 
plasmas (e.g.. iCioffi k, Joneslll980h . 

The intermediate situation, where two or three differ- 
ent RM components dominate within an individual ob- 
serving beam, is the most sensitive to the ambiguities 
discussed in this paper. The resulting complex interfer- 
ence patterns in A 2 space can give rise to erroneous RMs, 
and will increase the observed scatter preferentially on 
these angular scales. This situation will necessarily arise 
whenever the minimum angular scale of RM fluctuations 
is being approached. The only effective way to deal with 
this will be using Monte Carlo or numerical sim ulations 
(e.g., IGuidetti et ah! 120081 IGuidetti et all 120101 ) . where 
we expect there to be differences in shape between the 
input structure function and the observed structure func- 
tion on scales of the order of the beamsize. 

6. CONCLUSIONS 

We have presented our polarization analysis of compact 
radio sources observed with the WSRT at 350 MHz. Us- 
ing the observations of 585 sources in six fields, we com- 
puted a simple analytic model of the off-axis instrumental 
polarization (which can rise to several percent in q at the 
primary beam radius). After correction of the observa- 
tions using this model, only a small fraction of the sources 
were determined to have significant polarization at this 
frequency. By supplementing our observations with data 
from the NVSS, we have assessed the depolarization of 
our sample, finding the median depolarization ratio from 
1.4 GHz for the strongest sources to be P350/P1.4 < 0.2. 

We modeled the Faraday structure of seven sources us- 
ing various methods, including the traditional linear fit 
to x(A 2 ), as well as q, u vs. A 2 fitting to two simple de- 
polarization models - a foreground screen and two inter- 
fering RM components. In addition, we applied the novel 
Rotation Measure Synthesis and RM Clean techniques. 
A comparison of the RMs determined by various meth- 
ods has shown agreement in many sources, and yet fail- 
ure to reproduce the q, u observations casts doubt upon 
the validity of those solutions. In only one of the seven 
sources modeled, where depolarization from 1.4 GHz was 
not present, did the linear x(A 2 ) fit offer a solution that 
sufficiently reproduced the q, u observations. Of the re- 
maining six sources, RM Synthesis/Clean suggested mul- 
tiple significant (pm/poi > 0.5) RM components in three 
sources, while the two-component model found a signifi- 
cant secondary RM component in all six. Thus, a "char- 
acteristic" RM may be said to exist for any source, but 
the true Faraday structure may not always be adequately 
described by this alone. This point is well demonstrated 



by our detailed analysis of the southern lobe of the ra- 
dio galaxy 3C33. Previous studies, as well as our own 
linear x(-^ 2 ) fit and RM Synthesis/ Clean analyses have 
found a single, dominant RM of -7 rad/m 2 , in sharp con- 
trast to the q, u observations which strongly suggest two 
significant RM components near -3 and rad/m 2 . 

To further explore the possible shortcomings of the lin- 
ear x(A 2 ) fit and RM Synthesis/Clean methods, we have 
performed a few simple experiments. By constructing 
synthetic q, u spectra using the best fit two-component 
model for 3C33S, we find that RM Synthesis may place 
power at incorrect Faraday depths when multiple, closely 
spaced RM components interfere. In this case, both RM 
Synthesis and the linear x(A 2 ) fit find a consistent solu- 
tion, but one that does not agree with the known model 
inputs. The vulnerability of RM Synthesis is further 
demonstrated by a second experiment, which illustrates 
the role of phase in RM ambiguity. In this experiment, 
we show that two RM components, separated by more 
than the FWHM of the RMSF, may still yield an incor- 
rect solution under RM Synthesis depending upon the 
relative phase (i.e. intrinsic polarization angle) between 
the two components. A third experiment shows the dan- 
gers of a common assumption, that RM determinations 
made at high frequencies are sufficient. We show that 
the A 2 coverage must be broadened as much as possible 
to explore the true depolarization behavior of the source. 

With modeling of our WSRT observations and exper- 
iments on synthetic observations, we have touched upon 
some of the ambiguities that exist in rotation measure 
determinations. We caution that care must be taken 
when designing RM experiments and choosing one or 
more methods of analysis, stressing the importance of 
considering both degree and angle of polarization (or 
equivalently, q and u) over as wide a range of A 2 space as 
possible to see a more global picture of the polarization 
behavior and produce a more accurate description of the 
Faraday structure. 
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